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(Dated: February 6, 2008) 

We present a measurement of the cross section for Z production times the branching fraction to 
r leptons, (T-Br(Z — > t'^t~), in pp coUisions at -^s =1.96 TeV in the channel in which one r decays 
into ^v^iVt, and the other into hadrons + i^r or eVeVT. The data sample corresponds to an integrated 
luminosity of 226 pb~^ collected with the D0 detector at the Fermilab Tevatron collider. The final 
sample contains 2008 candidate events with an estimated background of 55%. From this we obtain 
(T-Br(Z T^T^) — 237 ± 15(stat)±18(sys)±15(lum) pb, in agreement with the standard model 
prediction. 

PACS numbers: 13.38.Dg,13.85.Qk,14.70.Hp 



Measurements of the Z boson production cross section 
times the leptonic branching fraction (a-Br) in pp colli- 
sions can be used to test standard model (SM) predic- 
tions. The (T-Br to e'^e~ and ii~ in pp colhsions has 



been measured by the UAl and UA2 collaborations at 
= 630 GeV ^l], by the CDF collaboration at Vs = 1.8 
TeV and ^ =1.96 TeV 0, and by the D0 collabora- 
tion at y/s = 1.8 TeV The Z boson branching ratio 
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to r+r~ has been measured with high precision by the 
CERN e~^e~ colhder (LEP) experiments These mea- 
surements are in good agreement with SM expectations 
and lepton universahty. We report here the first mea- 
surement of cr -Br (Z —> T+T~) in colhsions. This mea- 
surement provides a test of the SM as a deviation from 
the expected value would be an indication of anomalous 
production of r+r" pairs in pp collisions. It also veri- 
fies that the D0 detector can identify isolated t leptons 
in the energy range covered by Z boson decays, which 
could be critical in the search for non-SM signals such 
as supersymmetric (SUSY) particles in certain regions 
of the SUSY parameter space, or heavy resonances de- 
caying into fermion pairs with enhanced coupling to the 
third generation. 

The D0 Run II detector is fully described in jsj; a 
more succinct description of details relevant to this mea- 
surement can be found in 0. The Z — > t{~^ ^v^Vr)T 
candidate selection strategy focused on one r lepton de- 
caying to muon by triggering on the single muon using 
a three-level triggering system. The first level used the 
timing and position information in the muon scintillator 
system to find muon candidates. The second level used 
digital signal processors to form segments defined in the 
muon drift chambers. The third level used software algo- 
rithms executed on a computer farm to reconstruct tracks 
in the central tracking system and required at least one 
track with transverse momentum pT > 10 GeV. The in- 
tegrated luminosity of the selected sample is 226 pb~^ 
determined with a 6.5% uncertainty 0- 

After full reconstruction, the events were required to 
have an isolated muon with > 12 GeV and a r can- 
didate. The muon isolation required less than 4 GeV in 
the calorimeter in a cone R = \J (A0)^ -f (A?7)^ < 0.1 
(where is the azimuthal angle and 77 is the pseudora- 
pidity) around the muon, less than 4 GeV in an annulus 
0.1 < i? < 0.4, and fewer than three tracks (other than 
the muon) with pT > 0.25 GeV within R < 0.7. 

Most T leptons decay to one or three long lived charged 
particles plus up to three 7r° mesons that can be ob- 
served in the detector. The t candidates were found by 
constructing a calorimeter cluster made of all the towers 
with energy above a preset threshold around a seed tower 
within R < 0.5, keeping only clusters with > 5 GeV 
and E^"'''^ > 4 GeV, where E:^ {Eif'''') is the transverse 
energy with respect to the beam axis within R < 0.5 
{R < 0.3), and requiring rmsr < 0.25 (see Table cap- 
tion) and at least one associated track withpr > 1.5 GeV 
within R < 0.3. If there was more than one track, the one 
with highest pr was associated with the r candidate. A 
second track was added if the invariant mass calculated 
from the tracks was less than 1.1 GeV, and a third if the 
invariant mass was less than 1.7 GeV and the total charge 
was not ±3. Candidates with total charge zero were dis- 
carded. Finally, subclusters were constructed from the 
cells in the EM section of the calorimeter belonging to the 



TABLE I: Event pre-selection cuts 



Selection 



applied to the r-types 



only one fj, 
pt^ > 12 GeV 
/i isolation 

> 10(5) GeV 
^p^r'"'" > 7(5) GeV 
rmsr < 0.25" 
|<^^ -<^,i > 2.5 



all 
all 
all 

1 and 3 (2) 
1 and 3 (2) 
all 
all 

1 and 2 



Vms^ = \/Ei'=i[{^9^i)^ + (A-niy^]ET./ET, where i = 1, ..,n is 
the index of the calorimeter tower associated with the r-cluster; 
Arji and A(f>i are the rj and difference between the center of the 
T-cluster and calorimeter tower i. 

''Tt'l^f, = (E^ —E^^)/pif'', where E^^ is the energy deposited in 
a window of 5 X 5 towers (each tower of size <j) X r]=0.1 X 0. 1) around 
the T-track in the coarse hadronic (CH) section of calorimeter. 



T-cluster. The minimum Et required for an EM subclus- 
ter was 800 MeV. Three types of r candidates were iden- 
tified according to tracking and calorimetry information: 
1) single track with no subclusters in the electromagnetic 
(EM) section of the calorimeter (7r-like), 2) single track 
with EM subclusters (p-like), or 3) more than one asso- 
ciated track. No attempt was made to separate hadrons 
from electrons (which can contribute to both r-type 1 
and r-type 2). 

Additional requirements (which depend on the r-type) 
imposed on the selected events to enhance the signal-to- 
background ratio are shown in Tabic ^ The background 
increases rapidly with decreasing pt^ or decreasing E^. 
It is significantly lower for r-type 2 than for the other r- 
types, so a lower E^ cut is warranted for that r-type. The 
10^1 — 0t| > 2.5 cut takes advantage of the fact that most 
Z bosons have low px and thus the decay r leptons are 
back-to-back in (j). The longitudinal shape variable T^.J'^j, 
(defined in Table ^ caption) is used to remove misiden- 
tified muons because it has a distribution that peaks at 
much lower values for muons than for r leptons. 

The r leptons from a Z boson decaying to hadrons -I- 
i/r have average visible energy {E'^) of the order of 25 
GeV and need to be separated from a very large back- 
ground of jets. To further reduce the jet background, a 
neural network (NN) consisting of a single input layer 
containing several nodes (one for each input variable), a 
single hidden layer with the same number of nodes, and a 
single output node was used. A separate NN was trained 
for each type using a Monte Carlo (MC) sample of single 
r leptons uniformly distributed in Et and 77 and over- 
layed with a minimum bias event for signal |9|, and jets 
recoiling against non-isolated muons from data for back- 
ground. The NN input variables were chosen to minimize 
the dependence on the r energy and to exploit the nar- 
row width of the energy deposition in the calorimeter, the 
low track multiplicity, the low r mass, and the fact that 
r leptons from Z boson decays are well isolated. The NN 
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input variables were: 

1. profile ~ {Eti + Et.2 )/E^, where Et^ and are 
the Et of the two most energetic calorimeter tow- 
ers. Used for all r-types. 

2. caliso = {E:^ - Ei^°'''') / Ef" . A calorimeter isola- 
tion parameter used for all r-types. 

3. trkiso — Yjp^^^ /YipT^'''' , where p^^^ {Pt'''') the 
Pt of a track within a i? < 0.5 cone not associated 
(associated) with the r candidate. A track isolation 
parameter used for all r-types. 

4. {E^^'^^ +£;£;m2 ) /^r in a i? < 0.5 cone, where E^^'^^ 
and are the energies deposited in the first two 
layers of the EM calorimeter. A parameter used for 
r-type 1 to reject jets with one energetic charged 
track and soft tt" mesons. 

5. p^'"'*'^ /E^, where p'^''''^ is pt of the highest pr track 
associated with the r. Used for r-type 1 and 3. 

6. pTj;'^''^ /{EJp- caliso). A parameter used for r-type 2 
that measures the correlation between track and 
energy deposition in isolation annulus. 

7. ei2 = ^Wt*' ■ E^^' /E^, where E!^^ is the 
transverse energy deposited in the EM layers of the 
calorimeter. Used for r-types 2 and 3. 

8. 5a — ^ ( A(/)/ sin 0)^-1- ( Ar;)^ , where the differences 
are between Er-tracks and SEM-clusters. In the 
small angle approximation the observed r mass is 
given by ei2 • ETj, ■ 8a. Used for r-types 2 and 3. 

The dominant background is from multijet (QCD) pro- 
cesses, mainly from && events where the muon isolation 
requirement is met and a jet satisfies the r selection cri- 
teria. The other sources of background are W —* p^v 4- 
jets and Zj^* with one of the muons misiden- 

tified as a r lepton. The > 0.7 cut removed 70% of 
the M^/i^ background while keeping 98% of the expected 
Z/7* T^T^ events. The number of events that did not 
satisfy this criterion was used to estimate the background 
from misidentified muons remaining in the sample after 
the cut. 

The selected 29,021 events were separated into two 
samples: [i and r of opposite charge sign (OS), and /i 
and r of same charge sign (SS). The OS sample contains 
the signal. The SS sample is dominated by background 
and was used to predict the QCD background distribu- 
tions in the signal sample. From detailed studies of a 
sample of data with non-isolated muons, we established 
that this procedure is sound if one accounts for a small 
excess of OS over SS events that varies somewhat with 
the r-type. The correction factors (/i, where i denotes 
the r-type) were determined to be 1.06±0.06, 1.09±0.03, 



and 1.03 ± 0.02, by taking the ratio of OS to SS data in 
the non-isolated muon sample. There was no observable 
dependence of ji as function of Elj. , of NN output (NN) 
values, or of the muon parameters. A possible depen- 
dence of fi on the degree of isolation of muons coming 
from jets was checked by looking at the variation in the 
fi as function of muon px relative to the jet axis and 
varying the muon isolation. No variation was observed 
within the systematic uncertainties quoted. An overall 
3% systematic uncertainty was added for the extrapola- 
tion to the NN > 0.8 region. These factors do not fully 
account, however, for the contribution from W fiv + 
jets, which have a larger excess of OS over SS and dif- 
ferent distributions. The additional contribution of this 
channel to the signal sample is estimated from pythia 
MC samples. The MC is normalized using the OS 
and SS data with > 20 GeV, |(/)^ - (^r| < 2.0, and 

0. 3 < NN < 0.8 (in this region W ^ fiiy events dominate 
over the QCD background). The additional contribution 
to the background from W ^ tv events was ignored as it 
is a small fraction of the uncertainty on the contribution 
from W pv events. 

Figure n shows the NN distributions for each r-type 
(and the sum) for the signal sample, the predicted back- 
ground and the result of adding the predicted signal 
(from Z/7* ^ rr MC 0) to the background. Table |Il| 
shows the total number of events observed and predicted 
before and after the final cut NN >0.8. Distributions of 
background subtracted data are in very good agreement 
with those expected from Z ^ tt MC. Figure |21 com- 
pares the expected E^ and pt^ (adding all r-types) signal 
distributions to the predicted background distributions, 
and to the distributions obtained by subtracting the pre- 
dicted background from the signal sample distributions. 

The total event efficiency (ctot) summed over r-types 

1, 2, and 3 is 1.52% for M^r greater than 60 GeV . The 
total efficiency accounts for all losses due to branching 
ratios, geometrical acceptance, reconstruction and trig- 
ger efficiencies. It is corrected for the small difference 
between MC and data reconstruction cfRciencies. The 
contributions of the three r-types to the signal in the 
final data sample are 13%, 58%, and 29%. 

The cross section times branching ratio for Z/j* 
T+T- is given by A^signai/(eTOT • / Cdt) where iVsig„ai is 
given by the number of signal events and / Cdt is the 
integrated luminosity of the sample studied. iVsignai = 
865±55 (statistical uncertainty only) is the number of 
OS events of all r types after selecting the events with 
NN > 0.8, subtracting the estimated background (see 
Table ITl|. and subtracting the number of expected events 
in the sample with AIt-t less than 60 GeV (3.5%). 

The systematic uncertainties on the cross section mea- 
surement are listed in Table ITTTl The uncertainty (2.5%) 
due to the energy scale was estimated from the change 
in the acceptance when scaling the energy in MC events 
by the energy difference between MC and data (as de- 
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FIG. 1: NN output distributions for: (a) r-type 1, (b) r-type 
2, (c) r-type 3, and (d) the sum over all the r-types. 



TABLE II: Number of predicted and observed contributions 
to OS events by r-type before and after the A'^A'^ > 0.8 cut 





r-type 1 


r-type 2 


r-type 3 


Total 


QCD" 


1638 ± 107 


6001 ± 187 


6242 ± 153 


13881 ± 264 


Z/j* fifi 


33± 11 


67± 22 




100 ± 24 


W nu^ 


42 ±41 


151 ±93 


241 ± 114 


434 ± 153 




139 ±6 


700 ± 26 


335 ± 14 


1174 ±43 


Sum 


1852 ± 117 


7019 ± 214 


6818 ± 189 


15589 ± 309 


OS events 


1880 


6971 


7060 


15911 


AA > 0.8 


QCD 


196 ± 23 


280 ± 24 


508 ± 32 


984 ± 46 


Zf-y* fifi 


30± 10 


40± 13 




70± 16 


W —* fiU 


3±5 


17± 11 


38± 16 


58 ±20 


Z/'y* —* TT 


121 ± 6 


532 ± 21 


261 ± 11 


914 ± 24 


Sum 


350 ± 26 


869 ± 36 


807 ± 37 


2026 ± 57 


OS events 


355 


820 


833 


2008 



"The QCD background is estimated by multiplying the number 
of SS events by fi (described in the text). 

''The expected contribution is the number of events that must be 
added after subtracting the corrected number of SS events from OS 
events. 

■^The predicted number of Z/7* t~^t~ events is ba sed on a 
theoretical cross section of 257±9 pb for Mrr > 60 GeV E3 plus 
3.5% predicted from MC for the number of events expected with 
Mrr < 60 GeV. 



termined by the px imbalance in photon + jet events). 
The systematic uncertainty due to the NN performance 
(2.6%) was estimated by generating ensembles of Monte 
Carlo events in which the number of events in each bin of 
distributions of NN input variables was allowed to fluctu- 
ate by the uncertainties in the difference between MC dis- 
tributions and the background-subtracted data distribu- 



FIG. 2: [(a), (c)] and [(b), (d)] distributions after 

AA > 0.8 cut: (a), (b) estimated background (open trian- 
gles) and predicted Z ^ tt signal (histogram); (c), (d) back- 
ground subtracted data (open circles) and predicted Z —^ tt 
signal. 



TABLE III: Systematic uncertainties on cr-Br(Z/7* —> r^r ) 


Energy Scale 


2.5% 


NN 


2.6% 


QCD background 


3.5% 


Z/'y* nil background 


2% 


W —* 111/ background 


2.3% 


Zh* TT MC 


1.5% 


PDF" 


1.7% 


edata/eA/C 


2.1% 


Trigger 


3.5% 


Total 


7.5% 



"Efficiency uncertainty due to uncertainty in parton distribution 
function (PDF). 

^^data/^MC i^ the ratio of data to MC reconstruction efficiency. 



tions. The distributions of NN input variables arc in good 
agreement with those predicted adding Z/^* — > t^t~ 
MC and the estimated background; two are shown in 
Fig. El 

The QCD systematic uncertainty (3.5%) is due to the 
uncertainty in determining fi. The uncertainty in the 
Z/7* M"*"/^" and W ^ backgrounds (2.0% and 
2.3%) come from the statistical uncertainty in determin- 
ing their contribution, while the Z/^* t'^t~ MC sys- 
tematic uncertainty reflects limited signal MC statistics. 
The edata/cA/c is dominated by the uncertainty in esti- 
mating the difference in r-type 3 tracking efficiency be- 
tween MC and data using the the ratio of two- to three- 
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FIG. 3: Distributions for OS data, background and back- 
ground plus signal of two NN input variables before [(a), (b)] 
and after [(c), (d)]iViV > 0.8 cut: (a),(c) profile; (b),(d) cahso. 

prong events between background subtracted data and 
Z/'y* T^T~ MC in r-type 3 candidates. The uncer- 
tainty from differences in subcluster reconstruction, sin- 
gle isolated track reconstruction and muon isolation add 
up to about 1%. 

The trigger efficiencies were estimated using Z — > 
^JL^ data, the systematic uncertainty comes from the 
statistical uncertainty in that data; the uncertainties in- 
clude dependencies on rj and (f>. Systematic uncertainties 
from all other sources are less than 1%. Thus we obtain 

a ■ Br{Z/'j* tt) = 252 ± 16(stat) ± 19(sys) pb 

for Mrr greater than 60 GeV. The quoted statistical un- 
certainty is the uncertainty from OS and SS statistics 
(excluding the uncertainties on the correction factors). 
This yields, after removing the 7* contribution, 

cr • Br(Z tt) = 237 ± 15(stat) ± 18(sys) ± 15(lum) pb 

in good agreement with the NNLO standard model pre- 
diction of 242±9 pb 
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